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20.  Continued: 


the  target  and  interference  are  known. 

The  performance  of  this  matched  parameter  estimator  (MPE)  is  evaluated 
analytically  and  compared  to  the  conventional  estimator  and  to  MPE  simulation 
results  as  a  function  of  signal -to-noise,  interference- to-noise,  time  delay 
separation,  and  signal  and  interference  spectra.  In  addition,  a  potential  bias 
error  resulting  from  a  mismatched  MPE  is  also  evaluated. 
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INTRODUCTION 

The  conventional  approach  to  estimate  the  time  delay  between  two  sensors 
is  to  crosscorrelate  two  signals  and  search  for  the  peak  of  the  resultant 
crosscorrelation  function  [1,2].  However,  in  the  presence  of  interference, 
these  time  delay  estimates  are  biased. 

Two  potential  solutions  exist  for  removing  this  bias;  (1)  an  optimum 
multitarget  processor  [3],  which  requires  a  basic  reformulation  of  the 
estimator  and  (2)  a  matched  filter  estimator,  which  provides  additional 
multitarget  processing  capability  either  at  the  full-beam  or  split-beam 
level.  This  report  addresses  the  methodology  and  the  performance  predictions 
of  the  latter  approach  [4]  applied  at  the  generalized  crosscorrelator  (GCC) 
output. 


CONVENTIONAL  ESTIMATOR 

The  GCC  is  shown  ’n  the  left  portion  of  figure  1.  Assuming  that  the 
input  to  the  two  channels  is 

xi(t)  =  s(t)  I(t)  +  ni(t), 

and  ( 1 ) 

X2(t)  =  s(t+Di)  +  I(t+D2)  +  n2{t), 

where 

s(t)  =  primary  target  signal, 

I(t)  =  interfering  signal, 
ni(t)  =  channel  1  input  noise, 
n2(t)  =  channel  2  input  noise, 

Di  =  desired  target  time  delay,  and 

^2  =  desired  interfering  time  delay. 

Assuming  that  the  signal,  interference,  and  noise  are  joint  Gaussian, 
zero-mean,  uncorrelated  processes,  it  is  shown  in  appendix  A  that  the  GCC 
output  C(t}  for  large  averaging  time  T  may  be  represented  as 

2 

C(+)  =  Ri2(+)  ^  n{+)  =  X:  aj  Pj(+-Dj)  ^  n{t),  (2) 

j=l 

where  Ri2(+)  is  the  noise-free  crosscorrelation  function  between  the  two 
channels  and  n(+)  is  the  noise  component  whose  transform  has  zero-mean  and 
covariance 

N(f2)  «  Y  [G^lfi)  622^^1^  (3) 
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where  ai  and  32  are  the  respective  target  and  interference  power;  and 
P2  are  the  respective  normalized  target  and  interference  autocorrelation 
functions;  G^,  622*  and  G12  are  the  respective  auto-  and  cross-spectral 
densities  of  channels  1  and  2;  and  6(f)  is  the  familiar  impulse  function. 

The  conventional  estimator  simply  searches  for  the  global  peak  of  C(f). 
Denoting  and  as  those  values  of  ^  wnere  C{+)  and  Ri2{ peak, 
respectively,  it  is  shown  in  appendix  B  that  the  mean  and  variance  of  is 
given  by 


0 


0 


var(+^) 


where  is  the  solution  to 


00 

dw  w 

00 


G-, . 


12  +  w 

e  °  =  0, 


(4) 


(5) 


and 


w  =  2iTf , 


and 


Gij(f) 


.  Gi(f)  *  G2(f)  ♦ 


.  Gj(f)  ^Gjlf)  •G„^„^(f). 


-i2TrfD,  -i2irfDo 

=  G^(f)  e  ^  +  G^ff)  e 


where 

Gj  =  auto-spectral  density  of  the  primary  target, 

G2  =  auto-spectral  density  of  the  interference, 

Gnini  =  auto-spectral  density  of  the  channel  1  input  noise,  and 

*^n2n2  '  auto-spectral  density  of  the  channel  2  input  noise. 

When  interference  is  absent  (G2  =  0),  the  estimator  is  unbiased,  eg., 

+0  =•  ^0  =  and 
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var(D^) 


singl e  target 


2it 

=  r 


J  dw  ^^11^22  ~ 

^  00 


(6) 


This  represents  the  Cramer  Rao  lower  bound  for  estimating  the  time  delay  of  a 
single  target  in  uncorrelated  noise  and  will  be  used  later  for  normalization 
purposes. 

Figure  2  illustrates  the  typical  bias  characteristics  of  the  conventional 
estimator  in  the  presence  of  interference.  The  bias  error,  toBw*  is 
expressed  as  a  function  of  the  reciprocal  signal  bandwidth,  Bw,  and  is  a 
function  of  the  time-delay  separation,  (Di-02)Bw,  between  the  two 
interfering  targets  fcr  the  case  of  flat  spectra  having  unequal  bandwidths. 
Notice  the  large  bias  error  when  the  separation  is  small.  Conventional  bias 
reduction  schemes  tend  to  suppress  the  effects  of  the  interfering  source  by 
beam  shading.  We,  on  the  other  hand,  treat  the  time-delay  estimator  in  the 
presence  of  interference  as  a  two-parameter  estimation  problem  and  estimate 
rather  than  suppress  the  interference. 


MATCHED  PARAMETER  ESTIMATOR 


A  technique  for  removing  the  bias  error  is  to  match  the  assumed 
noise-free  GCC  output,  Ri2(f),  to  the  received  noisy  GCC  output,  C(+),  under 
a  least-mean-square  (LMS)  error  criterion.  Accordingly,  we  wish  to  minimize 
the  cost  function 


J(A^,  A^ ,  ^2^ 


r'1/2 

"2 

J  df 

SAj  ij(--yj)  -  c(.) 

"^1/2 

(7) 


where  T^  is  the  selected  observation  window  and  pj,  j  =  1,  2  are  our 
assumed  normalized  autocorrelation  functions  of  the  interfering  sources.  The 
matched  parameter  estimator  (MPE)  is  considered  to  be  matched  when  pj  =  pj 
and  mismatched  when  pj  ^  pj,  j  =  1,  2.  The  values  of  Aj,  yj,  j  =  1,  2 
at  the  minimum  value  of  J  represent  our  LMS  estimates  or  aj,  Oj,  j  =  1,  2, 
respectively. 

Because  Aj^  and  A2  may  be  determined  explicitly,  the  apparent 
four-dimensional  minimization  problem  suggested  above  may  be  reduced  to  a 
two-dimensional  maximization  scheme.  It  is  shown  in  appendix  C  that  J  is 
minimized  by  maximizing 

z(yi,y2)  =  (P1H2  *  ^2^1  -  2QHiH2)/F,  18) 


and  that  the  estimates  of  a^  and  32  are  given  by 


^1 


(P^  -  QH2)/F 


(9) 
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and 


where 


(10) 


r  1/2 

,.J 

■'l/2 


’1^^  -  yi)  ■  -^2^’ 


vJ 


1/2 


H..  =  I  dt  Pj(t  -  C(t),  j  =  1,  2, 

~^l/2 


/ 


1/2 


-2 


and 


Pj=J  d.  -p.  (t-y.)  .j  =  l,  2, 
~^l/2 


F  =  Pj  -  Q"  . 


(11) 


(12) 


(13) 

(14) 


Data  are  usually  received  at  only  discrete  values  of  t  rather  than  continuous 
values.  This  modification  is  easily  incorporated  in  the  above  expressions  by 
replacing  the  integral  over  t  in  (7),  and  (11)  through  (13)  by  a  summation 
over  the  discrete  values  of 


Expressions  (8)  through  (14)  represent  the  working  equations  of  the  MPE 
processor.  The  implicit  solution  of  (8)  for  the  estimates  of  and  O2 
requires  a  two-dimensional  peak-searching  algorithm  which  may  be  facilitated 
by  precomputing  and  storing  the  functions  P^,  P2,  and  Q.  Futhermore,  the 
process  of  locating  the  global  peak  of  2  may  be  accomplished  in  either  a 
tracking  or  acquisition  mode.  Without  apriori  knowledge  about  D]^  and  D2, 
one  needs  to  search  the  entire  (yi,  y2)  space  (acquisition  mode). 

However,  as  a  history  of  0]^  and  D2  is  established,  the  search  region  may 
be  reduced  and  (8)  may  be  solved  in  a  tracking  mode.  Once  the  global  peak  has 
been  found,  expressions  (9)  and  (10)  yield  estimates  of  the  primary  and 
interfering  target  powers,  respectively.  Finally,  the  residual 

Jm-in  =  ^2,  Oi,  52)  can  serve  as  a  goodness-of-f it  indicator. 

In  particular,  when  the  estimator  is  unbiased,  the  mean  value  of  is 

simply  the  total  noise  power  in  the  observation  window,  i.e.. 


mm 


di  n^(t) 


(15) 
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PERFORMANCE  PREDICTION 

It  is  shown  in  appendix  D  that  the  MPE  estimates  of  Di,  03,  ,  and 

33  are  unbiased  provided  the  GCC  averaging  time  T  is  large  and  the  MPE  is 
matched,  eg.,  pj  =  pj.  This  statement  applies  regardless  of  whether  the 
data  are  discrete  or  continuous  and  without  restriction  on  the  observation 
window  Tj . 

When  the  estimator  is  mismatched  (pj  ^  pj),  the  above  estimates  will 
be  biased.  When  the  mismatch  is  small  and  when  Ti->-<x.,  it  is  shown  (appendix 
0)  that  the  mean  bias  errors  of  Dj  and  O3  are  given  approximately  by 

b^U)  .aj'-Oj  Sj 

(16) 

b7(A)  =  62  -  S2. 

where 

$1  =  ^2  e  [^130^'''^^0  ■  %  *^20^  ~  ^220  ^10  ^loj^'^O 
^l[^120  ^20  *^20  ■  ^220^^'^'"o  ~  %  ^10^]  ^^0’ 

$2  =  i>2  e  [^130  Qio  ^10  "  ^llo(^'^''o  "  ^^O  ^2o)]/l'o 

^1  [^12o(^'^*'o  "  %  ‘lio)  ■  ^110  *120  '’2o]^’'o 

Weg  =  a^  E^  +  a2  ^2  e”’'*^^,  (18) 

2 

^0 "  ^110  ^220  ■  ^120 

^110  "  *^10  ■  ^10  ^10  ’ 

2 

^220  "  *^20  ■  *^20  II20  7*^0  ’ 

^120  "  ^120  %  I^IO  ll2o7*'o  ’ 
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K.q  =  /dt  ^  j  ^  2, 

*"0  "  ^10  ^20  "  ^0  ’ 

00 

PlO  -/_ff  il  . 

(19) 

P20  =£df  . 

-  OO  ^  ^ 

Qo  =£df  ^2  . 

Q^q  =  -i  J  df  =  -  Q2Q  , 

^120  ^1  ^2  ®  ’ 

W  =  2-rTf  , 

and  where  j&j,  oj  and  pj  are  Fourier  transform  pairs; 

Ej  =  0j  -  j  =  1,  2  is  the  spectral  mismatch  between  the  interfering 
targets;  and  a  =  Dj  -  D2  is  the  time  delay  separation  between  the  two 
targets. 

It  is  instructive  to  find  an  upper  bound  on  the  bias  error  (16).  Using 
(17)  and  (18),  one  may  show  in  a  strai ghtforward  manner  that 

00 

I  !^l!^  ~  ^220  ^0  1^2!^  “  ^110  ^0’  applying  the  Schwartz 

'  -.00 

inequality  to  (16)  and  using  the  above  result,  one  obtains 


The  variances  of  the  estimates  are  derived  (in  appendix  D)  under  the 
condition  that  the  MPE  is  matched  and  that  =  ».  The  results  are 
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var  (5^)  =-rVr/ 


var  (D^)  =  -5- 


"2  0 


and 


var 


var 


(i^)  =  \j  df 

— 00 

-  1  f°° 

(a^)  =  j  j  df 


00 

* 

?  /  * 

-iwDAS 

/df 

Si 

G22 

s 

"  (S2 

S  ^ 

VJ 

—  cc 

L 

00 

r 

-iwoA^I 

/df 

Gil  G22 

s 

^  (S2 

$2  e 

Vj 

—  00 

L 

- 

1  2 

( *  - 

-iwCy 

_^ll 

22 

"11 

+  I 

IS2  ''1  ® 

/J 

* 

n21 

1 2 

(  * 

-iwOi\^ 

Gii  G 

22 

''2I 

+ 

ri2  ''2  ® 

)\ 

> 

(21) 


(22) 


(23) 


(24) 


where 


Vj  ■  (P20  «J  -  Qo  *2  e'“)/Fo 


'2  -  'PlO  ‘2 


and 


■=11  Sin,- 


22  1  2 


r  *  -  ir  +  r 

^12  ”  ^  ^2  ^  * 


(25) 


A  worst  case  condition  occurs  when  =  D2  (a  =  0)  for  which  only  the 
spectral  characteristics  can  be  used  to  distinguish  between  the  two  interfering 
targets.  Generally  this  results  in  the  largest  variance  in  the  estimates. 

Under  this  condition,  Q^q  =  Q20  =  ^1  ^2  ^re  purely  imaginary  and 

VI  and  V2  are  real.  Expressions  (21)  through  (24)  reduce  to 


var  (Dj) 


A=0 


T"Tt/  ^^2  ^120  ■  ^1  ^220^^  ’ 

^1  ^0 


var  (D2) 


00 

2'  ^^1  ^120  ~  ^2  ^110^^  ’ 


A=0 


u  ^  T 
^2  *0 


(26) 


(27) 
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'  .=0  ° 


(2  8) 


and 


var  (a^j  ^/  df  (P  t.  ’  Qq  ^  ’ 

T  F  - 
a=0  ^ 


[29) 


where 


G=(G,*G.  )(G  +G  )^G  G 

1  ^'  ^2*^2 


ano 


=  2(G^  -  G^)"  ^  G_  . 


If  In  addition  to  a=0  the  interfering  spectra  are  identical  (^i  =  0^); 
then  Fq  and  \q  approach  zero  and  the  variances  (26)  through  (29)  approach 
infinity.  This  is  simply  an  indication  that  the  MPE  is  unable  to  distinguish 
between  the  two  interfering  targets  under  this  condition. 


RESULTS/DISCUSSION 

A  convenient  f i gure-of-meri t  for  comparing  the  performance  of  the 
conventional  ano  matched  estimators  is  the  normalized  degradation  ratio  (NOR), 
which  is  defined  as  the  total  rms  error  in  the  presence  of  interference 
dividea  by  the  root-mean-square  (rms)  error  of  the  conventional  estimator  in 
the  absence  of  interference.  The  latter  normalizing  term  is  simply  the  square 
root  of  expression  (6).  Note  that  NOR  >  1  and  that  large  values  suggest  that 
the  accuracy  of  the  primary  target  time-delay  estimates  is  seriously  degraded 
in  the  presence  of  an  interfering  target. 

Comparisons  of  the  NOR  for  conventional  and  matched  estimators  as  well  as 
the  corresponding  simulation  results  are  shown  in  figures  3  through  6  for  the 
case  of  flat  spectra  and  various  noise  powers.  In  each  case  the  MPE  is  assumed 
to  be  matched  and  the  signal  and  noise  bands  extend  from  0  to  B  (Hz/sec)  while 
the  interference  occupies  the  full  band  '  lures  3  and  4)  and  the  upper  half 
band  (figures  5  and  6).  Additionally,  qnal-to-interference  ratio  (SIR) 

is  1  dB  in  all  cases.  In  these  figurei  id  curves  represent  predicted 

results  for  the  MPE  as  obtained  from  (2  ;  while  the  triangles  represent 

corresponding  simulation  results.  Simi  he  dashed  curves  represent 

predicted  results  for  the  conventional  t  naior  as  obtained  from  (4),  (5), 
and  (6)  while  the  circles  represent  simuic  ion  results.  As  in  figure  2,  the 
time-delay  separation  is  expressed  in  terms  of  the  reciprocal  primary  signal 
bandwidth;  eg.  we  plot  NOR  as  a  function  of  (D^  -  D2)B.  Note  the 
excellent  agreement  between  the  predicted  and  simulated  results. 

When  the  spectra  are  identical  (figures  3  and  4),  the  standard  deviation 
of  the  MPE  estimates  approaches  infinity  at  zero  separation  because  the  target 
and  interference  are  indistinguishable  under  this  condition.  As  the 
separation  increases,  the  matched  NDR  decays  rapidly  to  a  fraction  of  the 
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conventional  NOR.  Note  the  separation  regions  where  the  conventional  NOR  is 
smaller  than  the  matched  NOR.  This  is  a  familiar  characteristic  that  often 
occurs  when  comparing  biased  (conventional)  with  unbiased  (MPE)  estimators. 

At  a  low  signal-to-noise  ratio  (SNR)  the  performance  improvement  is  not  as 
dramatic  because  the  random  error  is  much  larger  than  the  bias  error  and 
swamps  its  effect  (figure  4). 

Figures  5  and  6  show  an  even  greater  performance  improvement  because  the 
matched  estimator  is  able  to  resolve  the  target  and  interference  at  small 
separations  when  their  spectra  are  unequal.  As  before,  the  improvement 
degrades  at  low  SNR. 

A1 1  of  the  performance  improvements  in  figures  3  through  6  are  due  to  the 
fact  that  the  MPE  estimator  is  unbiased.  This  situation  occurs  only  when  the 
assumed  interfering  spectra  are  equal  to  the  actual  interfering  spectra.  The 
sensitivity  of  the  MPE  bias  error  due  to  a  spectral  mismatch  is  illustrated  in 
figure  7  for  the  case  where  the  primary  target  spectrum  is  known  (0-^B),  where 
the  actual  interfering  target  spectrum  (^2)  extends  from  B/2  to  B  Hz/ sec, 
and  the  assumed  interfering  target  spectrum  (^2)  has  the  same  bandwidth  but 
a  =^13  1/3^  center  frequency  mismatch.  This  figure  may  be  directly  compared  to 
figure  2,  which  shows  the  bias  error  of  the  conventional  estimator  for  the 
same  spectra  and  SIR. 

Examination  of  figure  7  reveals  that  the  bias  error  can  be  quite  large, 
but  is  smaller  when  the  center  frequency  of  ^2  up-shifted  rather  than 
aown-shif ted.  Comparison  of  figures  2  and  7  shows  that  the  MPE  peak  bias 
errors  are  2  to  3  times  smaller  than  those  of  the  conventional  estimator  and 
that  they  occur  at  a  smaller  time-delay  separation.  These,  as  well  as  other 
unpublished  results,  tend  to  suggest  that  peak  MPE  bias  errors  are  usually 
less  than  conventional  bias  errors  (even  for  complete  spectral  mismatch)  and 
that  the  smallest  peak  bias  errors  are  obtained  when  the  higher  spectral 
frequency  components  are  well  matched. 

CONCLUSION 

The  MPE  technique  yields  unbiased  time-delay  estimates  in  the  presence  of 
interference  provided  the  functional  form  of  the  interfering  spectra  are 
known.  In  comparison  to  conventional  estimators,  rms  errors  can  be  reduced 
several  orders  of  magnitude  particularly  at  high  SNR  and  INR.  When  the 
assumed  interfering  spectra  are  unequal  to  the  true  interfering  spectra,  the 
MPE  estimates  are  also  biased. 

In  practice,  apriori  information  about  the  primary  target  spectrum  must 
be  known  for  both  conventional  and  MPE  estimators.  Thus,  the  decision  to  use 
the  W^E  rests  on  what  apriori  or  measured  information  exists  about  the 
secondary  interfering  target  spectrum.  Simulation  results  tend  to  suggest 
that  peak  MPE  bias  errors  are  usually  less  than  corresponding  conventional 
bias  errors  particularly  when  the  higher  spectral  frequency  components  are 
well  matched. 
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ESTIMATOR 


Figure  1.  Conventional  and  Matched  Time-Delay  Estimators 
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Figure  3.  Oegration  Ratio  Versus  Time-Delay  Separation 
(SNR  =  0  dB,  INR  =  -1  dB;  S;0  ->  B,  I;B/2  B.  N:0  -»  8) 
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Figure  5,  Oegratlon  Ratio  Versus  Time-Delay  Separation 
(SNR  =  0  dB,  INR  =  -1  dB;  S:0  8.  I:B/2  B,  N:0  B) 
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Figure  6.  Degration  Ratio  Versus  Time-Delay  Separation 
(SNR  =  -10  dB,  INR  =  -11  dB;  S:0  B.  I:B/2  B,  N:0  B) 
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Figure  7.  MPE  Bias  Versus  Time-Delay  Separation 
(SNR  =  0  dB.  INR  =  -1  dB;  S:0  -»  B,  l:B/2  -*  B,  N:0  -»  B) 
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APPENDIX  A 

representation  of  the  generalized  crosscorrelator  (gcc)  output 


Referring  to  the  GCC  illustrated  in  figure  1  and  the  channel  inputs  (1) 
described  in  the  text,  the  GCC  output  V(+ )  is 


V( 


T/2 

=  / 


dt  Xj^(t)  X2(t~  ^) 


=  j 


(t)  X^(t-t)  dt. 


-T/2 


(A-1) 


where  II(y)  = 


1  It  I  <  T/2 


[0  otherwise 


Its  mean  value  is 

vni 


=  Y  /dt  IKy 


)  x^(t)  x^(t-t)'  = 


where  Ri?(  is  the  crosscorrelation  function  between  the  two  channels. 
The  covariance  of  v(  t)  is  defined  as 


(A-2) 


Cov(  :  v(  +^)  v(  t^)  -  v(  t^)  .  v(  t^)  • 


(A-3) 


Substituting  (A-1)  and  (A-2)  into  (A-3)  and  using  the  familiar  chain  rule  for 
joint  Gaussian  rv's,  one  obtains 


S'  ffi  "‘2  "(r)  "(r) 


[Rl^(t^  tj,)  *^2  '2  " 


+  -  ^2  ^  *2^  ^^12^^  "■  *^2  " 


(A-4) 


the  respective  auto  and  crosscorrel at  ion  functions  of  the  two  channels,  i.e., 

2 

R^j(n  =iqTtr77[t^  =  2  *  ^,n, (O’ 

j=i  ^  ‘ 


122(0  =  x^d)  X2(t-0  =  ^  aj  n.(0  -  R^  n  (O’ 

j=l  ^ 


(A-5) 


AO 
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‘)  =  Xj(t)  Xp(t-t)  =  ^  ‘ 

j=l 


where 

=  autocorrelation  function  of  the  target  (j  =  1,2), 
aj  =  RjiU)  =  power  of  the  target  (j  =  1,  2),  and 

Pj(!)  =  normalized  autocorrelation  function  of  the  target  (j  =  1,  2). 
Letting  v  =  t2  -  t2  in  the  integration  over  ti  and  noting  that 


dt^  1 1 


=  T  A(if), 


where  a(y)  = 


lul  <  T 


otherwise 


,  (A-4)  becomes 


Cov(‘p  ‘^)  =y  Jdu  :’.(y)  [R22(u)  R22(u  +  ^2  "  R22(^ 

(A-7) 

and  if  T  >>  ' +  |  ‘2! ,  then 

Cov(‘2,  ‘2)  '  T  -  ^2  ~  Rx2^'^^  ^2  “  (A-8) 

J  »  X) 

Now  we  wish  to  represent  the  random  GCC  output  v{^)  with  C(0  =  C(  ^  n(*' 

In  order  that  these  two  be  equivalent,  it  follows  that  CTH"  =  n{  ^ )  =  0, 

and  n(*^)  nl*^)  =  Cov(^^,  ^2)* 

Finally,  letting  N(f),  n(t);  Gii(f),  Rii(0;  G22(f),  R22(0; 

Gi2^f),  R12(0;  Gi(f),  Ri(0;  G2(f),  R2(0;  Gnimin,  RninilO; 
and  Gpppplf),  Rnon?(")  Fourier  transform  pairs;  it  follows  that 
N(f)  =^0'^and  the'^covariance 


=  /f 


-iZtrlf^  +  f^  t^) 


n(f^)  n(-2)- 


Substituting  (A-8)  into  (A-9)  and  performing  the  indicated  integrations,  it 
follows  that 


V'NTfp'  =  Y  [Gxx(fx)  ^22^^’  ^<^2  "  ^  ^^^2  ' 


{A-10) 


A-2 
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APPENDIX  B 

CONVENTIONAL  ESTIMATOR  STATISTICS 


The  conventional  ti;  ?  delay  estimator  locates  the  global  peak  of 
C(M  =  ^12(0  *  n(t)  or  the  zero  crossing  of  its  derivative.  Denoting  'g 
as  the  location  of  the  global  peak  of  Ri2(^)  and  expanding  C'(‘)  about 
by  means  of  a  Taylor  series  expansion,  one  obtains 

C’(t)  a  +  (t  -  f^)  C"(t^)  +  (B-1) 


where  the  primes  denote  differentiations  with  respect  to  t. 


Now  the  estimator  finds  the  value  of  for  which  C'(  +  ; 
this  value  as  our  estimate  ‘g,  we  obtain  from  (B-1) 


=  0.  Denoting 


C"(t  ) 


"  'o  -  " 


C'(t  ) 
0 


(B-2) 


r'l-o) 


where  the  Utter  approximation  is  valid  for  high  signal-to-noise  ratio. 
Expressing  and  n  in  terms  of  their  respective  Fourier  transforms  G12 
and  N,  we  obtain  (with  w  -  2TTf) 


r  iw+ 

C'(t„)  .  n'(*„)  - 

/df  (iw)  N  e  , 

•'—00 

(B-3) 

C’TV  =  '^"(^0^  = 

-  /df  w^  G^2  e  > 

mJ^CO 

(B-4) 

since  R‘(U)  =  0  (by  definition)  and  N  =  0.  It  follows  that  the  mean  of  the 
estimate  *0  =  'O  where  Iq  is  the  solution  to 


)  = 


k 


i  /df  w  G 


12 


=  0. 


(B-5) 


Using  (B-2)  and  (B-3),  the  variance  of  fg  is  given  by 


var(fo^  =  [C'(to)]^/[C"(t^)]^  , 


(B-6) 


where 


[C’(t„)]^  = 


■  -h 


df 


2  "'1  "2 


Nlf^) 


W2T 


(B-7) 


B-1 
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Substituting  expression  (3)  of  the  text  into  (B-7,'  anr)  pe^fc^ming  *he 
integration  with  respect  to  yields 


re,  ^|J‘  -I  jaf  [GjjG; 


Finally,  substituting  (B-8)  and  (B-4)  into  (B-6)  and  changing  the  integration 
variable  to  w  yields  expression  (4)  shown  in  the  text. 


B-2 
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Substituting  (C-7)  into  (C-2)  and  collecting  terms  yields 
J  =  Jq  +  [a2(P^  -  q2)  .  -  HjQ)  -  H^]/Pj. 


(C-8) 


Next,  differentiating  J  with  respect  to  A2,  one  obtains  aJ/aAp  =  0  when 
A2  satisfies 


A 


2  " 


-  H^Q 

“p7p~7  ’ 


(C-9) 


Substituting  (C-9)  into  (C-8)  and  collecting  terms  results  in 

J  =  Jq  -  z(y^, 


(C-10) 


where 


z(yi.  y2)  =  py 


(Hp  P,  -  H,Q)^ 

— TT  " 


(Pi  P2 


Q^) 


Pi  -  2Q  H2 

(Pj  P2  -  Q^) 


(C-11) 


(C-12) 


Since  q2  _<  Pj^  P2,  it  follows  that  Z  >  0;  and  since  Jq  >  0,  we  deduce 

that  J  is  a  minimum  when  Z  is  a  maximum.  Hence,  the  apparent  four-dimensional 

minimization  of  J(Ai,  A2,  yi,  yz)  can  be  achieved  by  maximizing  the 

two-dimensional  function,  z(yi,  y2);  and  A2,  Ai  can  be  obtained 

explicitly  from  (C-9)  and  ((i-7),  respectively. 

The  function  Z  has  an  interesting  geometrical  interpretation.  Dividing 
both  the  numerator  and  denominator  of  (C-12)  by  Pi  P2  yields 

2  2  2 

z  =  (n^  n2  -  2n^  n2  cos  e)/sin  ©,  (C-13) 


C-2 


A 
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2  2 

In  this  expression,  is  the  term  that  should  be  maximized  if  the  ref¬ 

erence  model  assumed  only  the  presence  of  the  primary  target.  This  is  evident 

2 

from  (C-8)  when  1^2  ~  Similarly,  r)2  should  be  maximized  if  one  assumes 

only  the  presence  of  the  interfering  source.  Considering  the  triangle  formed 
by  the  intersection  of  ni  and  n2  with  acute  angle  e;  note  that,  from 
(C-13),-^  sin  e  is  the  length  of  the  third  side.  It  follows  from  plane 
trigonometry  thatyfz' is  the  diameter  of  the  circle  that  circumscribes  this 
triangle.  Hence,  z  is  directly  proportional  to  the  area  of  the  c ircumscribed 
circle. 
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APPENDIX  D 

STATISTICS  OF  THE  MATCHED  PARAMETER  ESTIMATOR  (MPE)  ESTIMATES 


In  this  appendix,  we  derive  the  statistics  of  the  MPE  estimates  Di, 

D^,  ai,  and  32*  In  particular,  it  is  demonstrated  that,  without  any 
restriction  on  the  observation  window  Tj^,  the  mean  estimates  are  unbiased 
provided  the  MPE  is  matched,  eg.,  pj  =  aj.  Expressions  for  the 
approximate  bias  errors  of  and  62  are  then  obtained  for  large  Tj  when 
the  MPE  is  slightly  mismatched.  Finally,  expressions  for  the  variance  of  the 
estimates  are  derived  for  both  large  Tj  and  matched  conditions. 

The  bias  error  associated  with  estimating  the  coordinates  of  the  global 
peak  of  z(yi,  y2)  be  obtained  by  expanding  the  partial  derivatives  of 
z  about  the  point  (Di,  D2)  via  a  two-dimensional  Taylor  series.  Using  the 
last  subscript  0  to  denote  that  the  expression  be  evaluated  at  (Dj^,  O2) 
one  obtains 

zi{yi,  y2)  =  ^10  ^  "no  ^  '  ^2^  "120* 


where 


Z2(yi.  y2)  =  220  ^  "120  ^  (^2  -  ‘^2^  "220* 


Zj  =  32/3yj  ,  j  =  1,  2, 


«  J'  k  =  1,  2  . 

Now  the  search  algorithm  finds  those  values  of  yj  and  y2  for  which  the 
above  derivations  are  zero.  Denoting  these  values  as  our  estimates  of  D^ 
and  D2,  respectively,  and  solving  for  the  resultant  bias  error,  we  obtain 

"  2 
bi  =  Di  -  Di  =  {Zj20  ^20  “  "220  "10^ ^^"110  "220  "  "120^’  ' 


b,  =  D,  -  D,  =  (z 


120  "10  “  "no  "2o^^^"iio  "220 


"120^  ' 


(0-4) 


At  this  point,  let  us  separate  the  deterministic  and  random  parts  of  Hj 
and  z  by  defining 


h  .  (y  .  i 
J  -"j 


d+  Pj{ +  -  yj) 
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Pj( +  -  yj )  n( 0 


(D-6) 


such  that  Hj  =  hi  +  ej  for  j  =  1,  2.  Substituting  these  expressions 
into  equation  (8;  of  the  text,  we  obtain 

zly^.  y2)  =  6(^1.  y2)  y2^> 


where 


8 

=  (Pi  ^2 

"  P2  hJ  -  2  Q  hj  h2)/F 

(0-8) 

=  h^L  +  h 

=  (FM^  +  h^)/P^  =  (FL^  + 

(0-9) 

y 

=  2(Lej  + 

Me2)  P2  4  Pi  ■  n  ^2^/P 

(D-10) 

=  2(Le^  + 

Me2)  For  large  T, 

(0-11) 

L(y^, 

^2^ 

=  (P^  h^ 

-  Qh2)/F  =  (h^  -  QM)/Pp 

(0-12) 

M(y^, 

>2) 

=  (P^ 

-  Qh^)/F  =  (h^  -  QL)/P2, 

(0-13) 

F(yi. 

y2) 

=  Pi  P2  - 

(0-14) 

and  where  the  approximation  leading  from  (D~10)  to  (D-11)  is  valid  if  the  GCC 
averaging  time,  T,  is  large.  Note  that  hj  and  6  are  the  deterministic  parts 
of  Hj  and  z,  respectively;  and  that  y  ->  0  since  ej  :=  0. 


Now,  if  both  si gnal-to-noise  ratio  (SNR)  and  interference-to-noise  ratio 
(INR)  are  large,  we  may  replace  zjk  in  (D-3)  and  (0-4)  by  Bjk;  j,  k  =  1,  2. 
Substituting  zj  =  ej  +  j  =  1,  2,  the  bias  errors  reduce  to 


^1  *  ^®120^®20  ^20^  "  ®220  ®220  ~  ®120^ 
^2  *  1^®120^®10  ^10^  ■  ®110  ^®20  ®220  ~  ®120^ 


(0-15) 


I 


D-2 


I 
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where 

'j  ■ 

and  i 

=  3^6/ sy^  3y,^  ,  j,  k  =  1,  2  . 

Finally,  since  7  -  C>>  mean  values  of  these  bias  errors  are 

_  2 

bi  =  (0^20  ®20  ■  ®220  ®220  '  ®120^’ 

and  (D-16) 

_  2 

^2  ~  ^®120  ®10  ~  ®110  ®20^^^®110  ®220  '  ®120^ 

We  must  now  determine  the  various  derivations  of  B(yi,  yz)  evaluated 


D-3 
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and  di  and  are  the  signal  powers  of  the  primary  and  interfering 
targets,  respectively. 

Using  (0-12),  (0-13),  and  (D-17),  it  may  be  shown  that  for  the  matched  MPE 
Lq  =  "  ^2  ^ 


-^1  '^10 


■10  =  ~Wr 


'■20  "  '’26  "  ^^^20  ^*20^ 


^0  =  ^  P16  -  2Qio  Pio)  ) 


M2  = 


-32  F^o 

TFT 


Now  the  first  order  derivatives  of  B  may  be  obtained  from  (D-9)  and  are 
given  by 


=  L(2F^  +  LF^)/?^, 

=  L(2h^  -  2Q^M  -  LPp, 

82  =  M(2FM2  MF2)/P^. 

82  =  M{2h‘  -  20^1  -  MPp  . 


(D-20) 


(D-21) 


Substituting  (D-19)  into  the  above  expressions  at  yi  =  Di,  y2  =  D2,  we 
obtain  that  810  =  820  =  0*  It  follows  from  (D-16)  that  Bi  =  E2  =  0, 
eg.,  the  means  of  estimates  and  62  are  unbiased  when  the  MPE  is 
matched.  Also,  from  equations  (9)  and  (10)  of  the  text,  it  follows  that 
ai  =  Lq  =  ai  and  a2  =  Mq  =  32;  hence,  the  mean  power  estimates  are 
also  unbiased.  Notice  that  no  restriction  has  been  imposed  on  the  observation 
window,  T^,  and  that  the  MPE  remains  unbiased  for  discrete  data,  since  the 
integrals  are  then  simply  replaced  by  summations. 

The  second  order  derivatives  of  8  are  more  easily  obtained  from  (0-8). 
Although  tedious,  the  derivation  is  straightforward  and,  when 
yields 

2  2 

®110  "  ~^^1  ^110’  ®220  °  ~^®2  ^20’  ®120  "^^1  ^2 

where 

^10  "  *^10  ^'’16  %  ^10  ■  *^10  ^10  "  ^20  ’^16^'^'’^’ 


^20  *  "20  "  (P26  %  Q20  -  P20  Q20  '  PlO  ^20^/^)’ 

D-4 


pj  =  Pj’ 

Ii20« 


I 


1 

t 


1 
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When  the  MPE  is  mismatched  (dj  ^  pj),  all  the  above  estimates  will  be 
biased.  For  a  large  mismatch,  the  mean  estimates  of  Dj  and  D2  will  be 
given  by  the  simultaneous  solution  to  @1=6?=  0.  However;  when  the 
mismatch  is  small  we  may  use  (D-20)  and  (D-21)  to  derive  expressions  for  the 
approximate  bias  errors.  For  this  purpose  we  shall  assume  and  shift 

our  analysis  to  the  frequency  domain.  Using  Parseval’s  theorem  and  (0-5), 
(0-6),  (11),  and  (13),  we  may  express  hj,  ej,  Pj,  and  Q  as 


h  -  = 
J 

/df  6^2  ® 

y  J  = 

1,  2 

,  .  -iw(y .  -  D, ) 

ai  /df  ^ e 

+  32 

^  -iw(y.  -  Dp) 

e 

(D-24) 

e  .  = 

J 

Jdf  N  e 

,  j  = 

1,  2, 

(0-25) 

P.  = 
J 

/df  =  constant 

»  j  = 

1,  2, 

(D-26) 

Q  = 

jdf  e 

(0-27) 

where 


-iwD,  -iwDp 

Gj2(n  =  G^(f)  e  +  62(f)  e 

(0-28) 

-iwDi  -i202 

=  e  +  a2  ^2  ® 

and 

w  =  2itf , 


and  where  pj,  ij;  pj,  ^j  and  n,  N  are  Fourier  transform  pairs;  and  all 
integrations,  unless  otherwise  indicated,  are  integrated  from  -00  to  +=. 
Defining  Ej  =  ^ j  -  as  the  spectral  mismatch  between  the  actual  and 
assumed  normalized  auto-spectral  densities,  we  may  then  show  that  at  yj  = 


0-5 
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h 


10  =  "l  '’lO  "  ^2  %  '  ^1 

Jdf  w  weg 


*^l6  ■  ^2  ^10  ^ 


^20  ■  ^1  %  ^  ^2  '^20 


+  Jdf  (^2 


^20  ^  ^1  ^20  '^^2  ®'"*^  ’ 


TWA 


where 


weg  =  a^  Ej  +  32  E2  e 


-1  WA 


and 


A  =  Di  -  D^  =  time  delay  separation. 


and  where 

Qq  =  0^  ^2  ^ 


1  WA 


10 


■  k  *1 


«10  “  /“'  "  ‘1 


A  i-i”  ^  WA  n 

^2  ®  ■  '^20  ^20 


n20 


Jdfw‘‘  ^2  ® 


10 


/df 


-2Qo  Q^o  =  -^"20- 


Using  (D-12),  (D-13),  and  (D-29),  it  follows  that,  for  small  mismatch, 

Lq  =  a^  +  Jdf  weg  =  a^ 


Mo  =  a2  +  /df  weg  v,  a  ap  , 


where 


2  '  2 


''1  ■  ^^20  ^1  "  *^0  ^2  ® 


and 


V2  =  (Pio  62  -  Qo  Oi)/Fo  • 

Substituting  (0-29)  and  (0-31)  into  (0-20)  and  (0-21),  one  obtains 


(D-29) 


(0-30) 


(0-31) 


(0-32) 
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®10  “  wegliw^i^  +  Qjg  v^), 

®20  "  weg(iw^2  '^1  ^ 


(0-33) 


Because  the  mismatch  is  assumed  to  be  small,  the  second  order  derivatives  of  b 
may  be  replaced  by  their  matched  values  given  in  {D-22),  where  in  (0-23) 
reduces  to 


jn 

0 

1) 

'^10  ^10 

^20 " 

'^20  '  ^^20 

'^20^’^0’ 

(0-34) 

^20 " 

^120  ""  % 

*^10  ^20^''o  ’ 

since 

and  Kj  = 

|df  ,  j  =  1,  2  . 

Finally,  substituting  (0-22),  (0-33),  and  (D-34)  into  (0-16),  the  mean  bias 
errors  for  a  mismatched  MPE  are  given  by 

where 

2 

*0  "  ^10  ^220  "  ^120’ 

=  ^2  e  -  %  ^20*  “  ^220 

(D-36) 

^1^^120  ^20  ^^20  “  ^20 
S2  =  ^2  ®  [^120  ^10  ^10  ' 

(D-37) 

*  ^10^  -  ^110  ^^20  ^20^^''o  ’ 

weg  =  +  62  E2  , 
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-iwD, 


Jdf  N  e  ^  S, 


■  a,x 


ro 


/< 


df  N  e 


*  -iwDj 


~  a^x 


2*0 


(D-42) 


where  Vj  and  Sj  have  been  defined  previously.  It  follows  that 


var(O^)  =  b^, 


var  ( D 


"T? 


Ir^i  ^^2 

h  ^0 


)  e 


-12^D^(f^  +  f^) 


(D-43) 


Substituting  expression  (3)  of  the  text  into  (D-43)  and  integrating  once,  we 
obtain 


var(0,  )  =  y  K —  ic 
a2  x2  T  j 


df 


Gii  G22 


2  /  *  ■■’'"‘^1 
h\  Uhzh^ 


The  variance  of  D2  is  obtained  in  a  similar  manner  from  (D-42)  and  is 


var(D2)  = 


V-f 

xX  T  J 


L.  ^  T 

^2  0 


Gii  G22 


^2!^  +  IGi2  ^2  ^ 


-i  wOi 


Substituting  (D-25)  into  (D-39),  one  obtains 


var(d 


:)  =  /- 


*  -iwD, 


df  N  e 


/f  *  -iwD^  Y 
var(a2)  =  Mdf  N  e  V2  | 

where  and  V2  have  been  defined  previously. 


(D-44) 


(0-45) 


(D-46) 
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are  functions  of  a  so  that  all  the  statistics  are  functions  of  a.  Indeed,  a 
detailed  inspection  of  the  bias  and  variance  expressions  will  show  that  the 
mean  bias  errors  for  the  mismatched  MPE  are  odd  functions  of  a  and  the 
variance  expressions  are  even  functions  of  a. 
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